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We give a pedagogical review of a covariant and fully non-perturbative approach to study nonlinear 
perturbations in cosmology. In the first part, devoted to cosmological fluids, we define a nonlinear 
extension of the uniform-density curvature perturbation and derive its evolution equation. In the 
second part, we focus our attention on multiple scalar fields and present a nonlinear description in 
terms of adiabatic and entropy perturbations. In both cases, we show how the formalism presented 
here enables one to easily obtain equations up to second, third and higher orders. 

I. INTRODUCTION 

The relativistic theory of cosmological perturbations (see e.g. [1-7]) is an indispensable tool to interpret cosmological 
data such as the Cosmic Microwave Background (CMB) anisotropies, and thus to connect the scenarios of the early 
universe, such as inflation, to cosmological observations. Because the temperature anisotropies of the CMB are so 
small (5T/T ~ 10~ 5 ), considering only linear perturbations is an excellent first approximation and this is why most 
of the efforts devoted to the study of cosmological perturbations have dealt with the linear theory. 

However, with the rapidly increasing precision of cosmological data, we have now reached a stage where nonlinear 
features of primordial perturbations could be observationally accessible. In the last few years, this has motivated 
an intensive development of the theory of cosmological perturbations beyond linear order. Two main strategies have 
been considered. The first one is based on extending the traditional coordinate-based approach by adding second- 
order perturbations in the metric and in the matter fields and by computing the equations of motion for the metric 
and matter perturbations by brute force [6, 8-12]. The second strategy is based on a fully nonlinear treatment, but 
restricted to super-Hubble scales and working directly in a coordinate system [13-19]. This strategy is related to the 
so called separate universe picture that represents our universe, on scales larger than the Hubble radius, as juxtaposed 
Fricdmann-Lemaitre- Robertson- Walker (FLRW) universes with slightly different scale factors [20, 21]. 

The approach that we present here, based on [22-25], relies on a geometrical perspective. From a computational 
point of view, it can be seen as a "middle way" between the two main strategies discussed above. Inspired by the 
so-called covariant formalism [26, 27] for cosmological perturbations developed in [28-33] (see also the article by R. 
Maartens in this issue [34]), the fundamental idea is to use a nonlinear and covariant approach as long as possible and 
try to construct covariant quantities that mimic the traditional quantities that have been useful in the theory of linear 
perturbations (see also [35] for an alternative application of the covariant formalism to nonlinear perturbations). Since 
these quantities are tensors, they are intuitively much easier to understand from a geometrical point of view. Once 
these covariant objects have been defined and their evolution equation obtained in a fully nonlinear and covariant 
form, it is then possible to choose a coordinate system and to expand the quantities and equations to the desired 
order, in order to make concrete quantitative calculations. 

In this review we discuss two examples of this approach which are particularly illuminating. The first is the 
nonlinear extension of the familiar linear curvature perturbation on uniform-density hypersurfaces, usually denoted 
C, introduced in [36]. This quantity plays a crucial role in the linear theory, because it is conserved on large scales 
for adiabatic perturbations [37, 38], and therefore enables one to very easily relate cosmological perturbations in the 
very early universe and in the late universe, as long as the scales remain super-Hubble. The second example is useful 
in the context of multi-field inflation. In the linear theory, it is useful to decompose the perturbations of the scalar 
fields into (instantaneous) adiabatic and entropic modes, which are simply the projections of the perturbations along 
directions that are, respectively, tangential and orthogonal to the background trajectory in field space. Once again, 
it is possible to generalize this decomposition into adiabatic and entropic modes at the nonlinear level. This can be 
used to compute the evolution of nonlinear perturbations and the resulting non-Gaussianities. 

This review is divided into two parts. In the next section we consider the treatment of a cosmological fluid, while 
the third, and final, section is devoted to the case of several scalar fields. 
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II. COSMOLOGICAL FLUIDS 
A. Covariant approach 

In this section we consider a single perfect fluid, characterized by a comoving four-velocity u a (u a u a = —1), a 
proper energy density p and a pressure P. The energy-momentum tensor associated to the perfect fluid is given by 

n = (p + P)u a u b + Pg a b . (1) 

To fully characterize the fluid, one needs an equation of state relating P to p and, possibly, to other physical quantities 
if the fluid is not barotropic. 

The spatial projection tensor orthogonal to the fluid velocity u a is defined by 

hab =9ab+U a U bl (h a h h b r = h a cl h a b U b = 0). (2) 

It is also useful to introduce the familiar decomposition 

V 6 u a = cr ab + u ah + ^Qh ab - a a u b , (3) 

where one finds on the right-hand side the (symmetric) shear tensor <j a b 1 the (antisymmetric) vorticity tensor uj abl 
the volume expansion 8 = V a w a and the acceleration a a = u c \7 c u a . 

The integration of along the fluid world lines with respect to the associated proper time r, 



a = iy dre , (4) 



can be used to define, for each observer comoving with the fluid, a local scale factor S — e a . It follows that = 3d, 
where the dot of a scalar quantity denotes its derivative along u a , i.e. a = u a V a a. This quantity a, which can 
be seen as a (covariantly defined) local number of e-folds, plays a crucial role in the definition of our covariant and 
nonlinear extension of the familiar linear curvature perturbation on uniform-density hypersurfaces, usually denoted 
£, introduced in [36]. 

Our starting point is the conservation of the energy-momentum tensor, 

v„n = o, (5) 

which yields, after substituting (1) and projecting along u b , the continuity equation 

p + 3a(p + P) = 0. (6) 

The spacetime gradient of this expression can be written as 

(d aP y + 3 (d a a)- (p + P) + 3ad a (p + P) = 0, (7) 

where the dot acting on a spacetime gradient, which is a covector, is defined as the Lie derivative (see e.g. [39]) with 
respect to u a . For any covector X a , this means 

X a = C u X a = u b d b X a + X b d a u b . (8) 

More generally, a dot acting on any tensor will denote the Lie derivative of this tensor with respect to u a . This 
definition is compatible with our notation for the scalars since the Lie derivative coincides with the derivative along 
u a for scalars (and we have used the property d a (p) — (d a p)' to obtain (7)). 

After some simple manipulations, one finds that Eq. (7) is equivalent to [22, 23] 

(a = -177-^(da P --d a p) , (9) 



3(P + P) \ P 

where the left-hand side consists of the time derivative, i.e. the Lie derivative with respect to u a , of the covector 

( a = d a a - %d a p. (10) 
P 
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This covector can also be re-expressed, using (6), in the form 



(a = d a a + 

and, if w = P/p is constant, this is a total gradient since 

1 



Ca = d a 

On the right-hand side, the quantity 



3(1 + w) 



In p 



d a P 



Z(P + P) 



(w = P/p = const). 



(11) 



(12) 



r a = d a p - -d aP 
p 



(13) 



is a nonlinear generalization of the non-adiabatic pressure. It vanishes for purely adiabatic perturbations, for instance 
when the pressure P is solely a function of the density p. 

Eq. (9) has a form very similar to the conservation equation for £ of the linear theory, which will be rederived in 
the next section. At the non-linear level, our covector ( a is closely related to the non-linear perturbation introduced 
in [19], 



3Jp p + P 



(14) 



where p(t) is the homogeneous background density, while p(t,x) is the local energy density and N — f dt<d/3 is the 
number of c-folds, defined with respect to the time coordinate t. A similar non-linear quantity has been introduced 
earlier in [17]. The advantage of the definition (10) is that one does not need to introduce a coordinate system or to 
restrict oneself to super-Hubble scales, and that its evolution equation (9) is exact, fully non-perturbative and valid 
at all scales. Moreover, it is worth stressing that this equation is a direct consequence of the conservation of the 
energy-momentum tensor and is independent of the underlying theory of gravitation. 



B. Link with the coordinate approach 

We now relate the covariant approach with the more familiar coordinate based formalism [1—7]. We first examine 
the linear perturbations and consider later the perturbations at second and third orders. 



1. Linear theory 

A cosmological spacetime closed to FLRW geometry can be described by the perturbed metric 

ds 2 = -(1 + 2A)dt 2 + 2aB,dx l dt + a 2 ( 7ij + H^dx'dx 3 , (15) 

where a — e a is the background scale factor. As usual, one can decompose the linear perturbations into so-called 
scalar, vector and tensor modes, 

Hij = - 2ip 7ij + 2ViVjE + 2V (< £$ + 2Ej., 

where BY and EY are transverse, i.e., ViB Vl = = \7iE Vl , and Ej- is transverse and traceless, i.e., ViE Tl i = 

and -) l ^Ej- = 0. Here V» denotes the three-dimensional covariant derivative with respect to the homogeneous spatial 
metric 7y (which is also used to lower or raise the spatial indices). 

The corresponding matter content is a perfect fluid with perturbed energy density and pressure, p(t,x l ) = p(t) + 
5p(t,x l ), P{t,x l ) = P(t)+5P(t,x l ), and four-velocity 

u" = G" + <Su", 6u>" = {-A,v i /a}, v k = V k v + vY, (17) 

where vY is transverse, ViV Vl = 0. 
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It is easy to check that the covcctor ( a automatically vanishes in the background. At linear order, its spatial 
components are simply [22, 23] 

Ci = diC, C = 6a- j5 P , (18) 

where H denotes the Hubble rate, H — a' /a — a' and we recall that a prime denotes the partial derivative with 
respect to the cosmic time t. The quantity ( defined above differs from the familiar definition [36] 

Cb = -^-H% (19) 
P 



1 2n . 1 f ,. V 2 v 



because 5 a differs from —ip. They are related by [23] 

5a = -i> + ry' Z E+ : - I dt— , (20) 
o o J a 

and they coincide when the gradient terms are negligible. 

The components of the non-adiabatic term T a = d a P — (P/p)d a p can be deduced directly from the components of 
Ca by substituting P to a. One finds 

= diT, T = 5P- j5p. (21) 



Putting all together Eq. (9) gives, at first order, 



«--& + iV(J? + ./.)-^(*p-^). 



(22) 



For adiabatic perturbations, the right-hand side vanishes and the quantity £ is conserved at all scales because it comes 
directly from the covector £ a . By contrast, (b is conserved only when the gradient terms are negligible. 



2. Second-order perturbations and beyond 

The covariant approach turns out very useful to go beyond linear order in perturbation theory. Any function X 
can be decomposed in the form 

X(t,x l ) = X(t) + 5X {1 \t,x l ) + 5X {2) {t 1 x t ) + 5X (3 \t 1 x t ) + ... , (23) 

where a quantity with a superscript is the n-th order contribution to X in the perturbation theory expansion [9] 
(note that here we do not follow the convention of including a numerical factor 1/n! in front of the n-th contribution). 

Expanding Q — d- L a — (a/p)dip at second order one finds that, in contrast with the first order expression, £( 2 ) 
cannot be written as the gradient of a scalar perturbation. Indeed, after some manipulations one finds, at second 
order [22, 23], 

C f) = a,c (2) + 5 fd^\ (24) 

with 

C (2) - ^ (2) - f V 2) - j^'5p + jL Sp 5p> + J_ V. (25) 

On large scales, using 

a « In a - tp - i/> 2 , (26) 

one can show that £( 2 ) is directly related to the conserved second-order quantity defined by Malik and Wands [40]: 

~ Cmw ~ (Cmw) 2 ( see a l so t ne discussion in [41]). We use the symbol w to denote an equality only valid on large 
scales. 
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Now that we have identified a second-order expression for £, we can expand Eq. (9) at second order. On making 
use of the second line of Eq. (8) to reexpress the Lie derivative along u a in terms of the coordinate time derivative, 
and retaining only terms at second order, we have 



(i) diV J 



(27) 



Finally, on making use of Eqs. (18), (24), and (25), and that 9 = 3(1 — A)a' up to first order, we can explicitly write 
the conservation equation (9) up to second order and on all scales: 



A2), = H _ r (2) _ 1 _ r (i) c (iy _ ?L d .Ai) 
S p + P p + P S a * 



(28) 



The definition of can be read from the second-order expression of by substituting P to a. For adiabatic 
perturbations, we find that at second order the scalar variable £ is conserved only on large scales, when the last term 
on the right-hand side of Eq. (28) can be neglected. 

One can extend straightforwardly the procedure presented here to higher orders in the perturbation expansion 
[42, 43]. Expanding Eq. (10) at third order, after some manipulations one can write 



c (3) = *c (8) + |^c (2) ' + -^-d^' 



6pW 



Sp^_ 
2p>'< 



■*c (1)ft 



(29) 



with 



i 

p 



H 

7 



2p' 2 



bp 



1 (H 



(30) 



Sp\ 



Its evolution equation up to third order can be found by expanding Eq. (9). On large scales it reads [43] 



c (3) ' ~ -- 



H 



p + P 



,(3) 



1 _ rr (1 V (2)/ + rMfW) - 1 

p + P [ + ^ ' 3H(p + P) 



/\2 



r(C') 



(31) 



where T^ 3 ) is defined in the same way as C' 3 ^ in Eqs. (29) and (30) with the replacements a — >• P, ( — > T. 

In [42] one can find a discussion on perturbations beyond the third order, as well as a general expression for a 
quantity (( n \ defined at any order n, which is conserved on large scales for adiabatic perturbations. 



3. The issue of gauge-invariance 

One can check that, on large scales, £^ behaves as a gauge-invariant quantity. Given a vector £ a = J2 n ^ the 
gauge transformation it generates is defined by the transformation law of tensors [6] (whose coordinate functions are 
a particular case), 

T^e £ «T. (32) 

Since Ca vanishes at zeroth order, ( a is automatically gauge-invariant at first order, according to the first expression 
above. However, Q is not gauge-invariant at second order and the corresponding gauge transformation is given by 

Cf^CP+^CP-Cf+^oCP, (33) 

where for the last equality we have neglected the terms of higher order in spatial gradients, concentrating on large 
scales. By noting that at first order Sp/p' -» Sp/p' + it is easy to see using Eq. (24) that d^ 2 \ or cquivalently 

C^ 2 -*, is gauge-invariant at second order, on large scales. This can also be checked by directly applying the gauge 
transformation (32) on (( 2 ) defined in Eq. (25). This analysis can be repeated at third [43] or any higher order n [42] 
and one can show that the quantity is gauge invariant on large scales. 
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C. Dissipative and interacting fluids 

It is possible to extend our covariant approach to the case of dissipative and interacting fluids, as discussed in detail 
in [24]. In that case, the energy- momentum tensor is more complicated, 

Tab = pu a u b + Ph ab + q a u b + q b u a + Tr ab , (34) 

where the energy flow q a and anisotropic stress it ab satisfy 

q a U a = 0, TT ab = TT ba , TT a a = 0, TT ab U b = 0. (35) 

The continuity equation in this case reads 

p + e(p + P)=V, V = -(7r ab a ab + D a q a + 2q a a a ), (36) 

where D a denotes the spatially projected covariant derivative, orthogonal to the four- velocity u a . For a generic tensor 
T, the definition is (see e.g. [28]) 

D a T b . c ;-=h a d h b e ...hf...V d Tjr. (37) 

As a consequence of the presence of the dissipative terms on the right-hand side of (36), the evolution equation for £ a 
is slightly more involved and now reads 

3rV 2 

(a = — (r a + S a ) , (38) 

p 

where 

^ a = D a p-^D a p, = (39) 

This means that dissipative terms affect the evolution of the covector ( a in addition to the non-adiabatic term. 



III. COSMOLOGICAL SCALAR FIELDS 

We now consider the situation where matter is composed of several scalar fields, concentrating on the case [25]. We 
study N scalar fields minimally coupled to gravity with Lagrangian density 

£=-^V-%if), (40) 

where V is the potential and the summation over the field indices (/, J, . . .) is implicit. For simplicity we assume here 
canonical kinetic terms but the nonlinear extension to a large class of models with non-canonical kinetic terms, such 
as those studied at the linear level in [44], has also been developed in [45]. The energy-momentum tensor derived 
from the above Lagrangian reads 

T ab = da^idbV 1 - \g ab (^/dV + 2V) . (41) 

Given an arbitrary unit timelike vector field u a , it is always possible to decompose this energy momentum tensor in 
the form (34) with, in our case, 

P = ^(W J + V^V) +V, P=\(vi<P I -\D a <piD a <p^j -V, (42) 

q a = -tpiD a (p I , TT ab = DaLpjDbif 1 - ^h ab D c (p I D c (p I . (43) 



The evolution equations for the scalar fields are the ./V Klein- Gordon equations given by 

dV 
dtpi 



-V a VV/ + — = Cpi + Qipi + V VI - D a D a <p! - a a D aVl = 0, (44) 



where the second expression is obtained by using a decomposition into (covariant) time-like and space-like gradients 
defined with respect to u a . 
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A. Adiabatic and entropy covectors 



For simplicity we now restrict our analysis to the case of two scalar fields, which we will denote by <f> and \- In the 
two-field case it is possible to introduce a particular basis in the field space in which various field dependent quantities 
are decomposed into so-called adiabatic and entropy components. In the linear theory this decomposition was first 
introduced in [46] for two fields. This decomposition is particularly convenient to follow the time evolution of the 
curvature perturbation, which is sourced by the entropic perturbations. For the multi-field case it is discussed in [47] 
in the linear theory and in [48] in the nonlinear context. 

In our case, the corresponding basis consists, in the two-dimensional field space, of a unit vector e\ defined in the 
direction of the velocity of the two fields, and thus tangent to the trajectory in field space, and of a unit vector e\ 
defined along the direction orthogonal to it, namely 

e i = (4>/<j,x/<j) = (cos 0, sin 0) , e[ = (-\/ = (-sin0,cos0), (45) 

where we have introduced the formal notation a = (<ft 2 + x 2 ) 1//2 . Note that this notation can be misleading as, in 
general, in the nonlinear context a is not the derivative along u a of a scalar field a. 

The angle 0, which characterizes the rotation angle between the initial field basis and the adiabatic/entropic basis, 
is here an inhomogeneous quantity which depends on time and space. By taking the time derivative of the basis 
vectors e£ and ef , we get 



'-eel, el = ~K- (46) 



e 



Making use of the basis (45), one can then introduce two linear combinations of the scalar field gradients and thus 
define two covectors, respectively denoted by a a and s a , as [25] 

a a =eld a fi = cosd d a (j) + s'md d a x, (47) 
s a =eld a (fi = -sm6d a cf) + cosed a x- (48) 

We will call these two covectors the adiabatic and entropy covectors, respectively, by analogy with the similar definitions 
in the linear context [46]. Whereas the entropy covector s a is orthogonal to the four-velocity u a , i.e., u a s a = 0, this is 
not the case for a a which contains a "longitudinal" component: u a a a = &. It is also useful to introduce the spatially 
projected version of (47) and (48), 

0~a = e i D afl =0~a+ ™ a , S 1 EE e^D^j = S a . (49) 

The "adiabatic" combination of the Klein-Gordon equations, i.e., the contraction of (44) by e£, yields 

a + Ga + Vv = V a < -Y (sh (50) 

where we have defined 

V a = eiV VI , Y{s) = ^(s a + eai)s a . (51) 

Since e\(pi = 96-, the "entropic" combination of the Klein-Gordon equations, i.e., the contraction of (44) with el, 
gives simply 

&9 + V s = V a s a +Y (a) , (52) 

with 

K s = elV Vl , Y (a) ee I(s„ + . (53) 

In summary, the Klein-Gordon equations for the fields <p and x have been replaced by the equations (50) and (52), 
whose left-hand side has exactly the same form as the homogeneous equations in a Friedmann universe. However, 
these equations capture the fully nonlinear dynamics of the scalar fields and their right-hand side contains nonlinear 
(quadratic) terms represented by Yr s \ and Y( CT ) , sourcing the adiabatic and entropy equations respectively. 

We now consider the evolution equations for the covectors a a and s a , or rather for their spatial projections. The 
adiabatic evolution equation is given by 

(aa) 1 - + Qi&a) 1 + aD a Q + (v^ + e^j a 1 - D a (VV-) =(e-^j s a 

+ (0 - V as + 60) s a - D a Y {s) , (54) 
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where we have introduced the notation 

Vera — e a e cjV,(pi<pj 7 V,ss = G, s e s V,tpi <pj j V.sa = e s e iT^,'PI<PJ 7 (55) 

for the second derivatives of the potential. Note that the spatial projection and the time derivative do not commute 
for <7 a . The entropic evolution equation reads 



s a - * (* + V a )s a + (V ss - 9 2 )s a - D a (V c s c ) = -IBi&a) 1 - 
a 



-{v + V a )-6-V c 



a 



°i + D a Y (a) , (56) 



where the covectors s a and s a are purely spatial, i.e. (sa) 1 ' — s a and (sa) 1 - = s a . 

Starting from the fully nonlinear Klein-Gordon equations, we have thus managed to obtain a system of two cou- 
pled equations (54) and (56), which govern the evolution of our nonlinear adiabatic and entropy components [25]. 
Remarkably, they are rather simple and they look very similar to the linear equations for the perturbations 5a and 5s 
(see [46]). Furthermore, since these equations are exact, it is straightforward to expand them up to second or higher 
orders, as we will show explicitly in Sec. IIIC. 

B. Generalized curvature perturbations and large-scale evolution 

Let us now turn to the generalized curvature perturbation ( a , which we introduced in Eq (10) for a cosmological 
fluid. Although a system of several scalar fields cannot be described as a perfect fluid, its energy-momentum tensor 
is of the form (34) with 



p=^ + n) + K P = I(^-ln)-K 



(57) 



where we have defined 



q a = -aa^ , ir ab = Il ab - \h ab Jl, (58) 



Il ab = a^a^ + s a s b , I! = <t^<t ±c + s c s c . (59) 



The quantity £ a for the two-fluid system therefore satisfies an evolution equation of the form (38). 

In the linear theory, an alternative quantity to describe the primordial perturbation is the comoving curvature 
perturbation. A natural nonlinear extension of this quantity in our formalism is (see also [49] for an alternative 
definition) 

TZ a = -D a a - — — rj-qa . (60) 

As in the linear case, this expression can be directly related to the analog of the Mukhanov- Sasaki variables [50, 51], 
if we generalize the latter by defining 

Q{ ee £>V - ^-D a a (61) 

for each scalar field. The comoving curvature perturbation lZ a is then proportional to the adiabatic combination of 
the Qi, i.e. 

Qa = e aI Q{ -D a a, (62) 

a 

since one can write lZ a = a LpiQ I a /(cpj(p J ). 

Comparing the definitions of ( a and lZ a shows that they satisfy the relation 

Ca+^a = -"e a , e a = D a p - -^er^, (63) 
P v 

where the covector e a can be interpreted as a covariant generalization of the comoving energy density perturbation. 
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Let us now concentrate on the large scale limit by neglecting terms of higher order in spatial gradients. Using 
the energy constraint u a G ao u h = 8nGp and the momentum constraint u b G 0C h c a = 8nGq a , it can be shown that the 
comoving energy density perturbation, e a , defined in Eq. (63), can be neglected on large scales if the shear can also be 
neglected in this limit. This is the case in an expanding perturbed FLRW universe, where the shear rapidly decreases 
on large scales. 

Then, from their definition (49), a 1 and s a are first-order quantities with respect to spatial gradients while the 
scalars Yj- S ) and Y^ are second order, since they are quadratic in and s a (or their time derivatives). Hence, the 
right-hand side of Eq. (50) and of Eq. (52) can be neglected on large scales. Taking into account these simplifications 
and neglecting on large scales e a and the terms containing II, one finds that the evolution equations in Eqs. (54) and 
(56) simplify to give 

{da) 1 - + OK) 1 + &D a Q + {V aa - e 2 )^ « 2(0s a y - 26^s a , (64) 

a 

and 

Ma + 6s a + (V ss + 38 2 )s a w0. (65) 

While the entropy mode appears as a source term in the adiabatic equation, the entropy equation involves only the 
covector s a . 

On large scales, neglecting e a in Eq. (63) implies that the uniform density and comoving perturbations ( a and lZ a 
coincide (up to a sign), 

(a+Ka^ 0. (66) 

Moreover, in the same limit, the source term in the evolution equation for ( a , Eq. (38), is dominated by the gradient 
of the potential V along the entropic direction, so that 

2 9 2 

(a « -z — V>„. (67) 

3 p 

This equation is sufficient to describe the evolution of the adiabatic mode in the large scale limit, making Eq. (64) 
redundant. The analogs of the spatial components of Eqs. (67) and (65) have also been derived in [48]. 



C. Link with the coordinate approach 

1. Linear theory 

We now relate our covariant approach with the more familiar coordinate based formalism, starting with the linear 
theory. We use the metric and notation of Sec. II B 1 and the scalar fields are decomposed into a background and a 
perturbed part, ip!(t,x l ) — <fii(t) + S(pi(t 7 x l ). We now need to specify the components of the unit vector u a , which 
defines the time direction in our covariant approach. In contrast with the perfect fluid case, where u a has a physical 
meaning, the unit vector u a is here arbitrary. In a given coordinate system it is convenient to choose w M such that 
Ui = at first order. This implies that the components of w M and of the "acceleration" vector are given respectively 

by 

u M = {1 — A, -B i /a}, a"={0,V'i/ fl 2 }. (68) 

The background equations of motion can be deduced immediately from Eqs. (50) and (52) and read 

a" + 3Ha' + V a = 0, a'0' + V s = 0. (69) 

From their definition, one finds that the spatial components of a a and s a , at linear order, can be expressed as 

Sen = diScr, (5cr = cos 6 5(f> + sm8 S\, (70) 
Ssi = diSs, Ss = cosOSx — sm8S(j), (71) 

which coincide with the notations of [46]. 
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Linearizing the evolution equations for a a and s a , eqs. (54) and (56), one easily obtains the linearized equations for 
8 a and 8s. The latter can be written as 

6s" + 3H8s' + (V ss + W' 2 )8s - -^V 2 8s = -2^-5e , (72) 

a z a' 

where on the right-hand side appears the first-order comoving energy density perturbation 8e, defined by 

8ti = di5e, 8e = 6p — rrjSa, (73) 

which follows from the definition (63) of e a . The quantity 8e can in general be neglected on large scales, since it is 
second order in spatial gradients, as can be seen from the relativistic Poisson-like equation. 

Let us now turn to the adiabatic equation. Since 8a is not gauge-invariant, in contrast with 8s, it is more useful to 
consider the gauge invariant Sasaki-Mukhanov variable Qsm, defined as [50, 51] 

-i 

Q SM = Sa + —V- (74) 

Note that the above traditional definition coincides with the scalar quantity Q that can be extracted from our definition 
of Q a given earlier in Eq. (62), via Qi = diQ, only in the large scale limit, because ip and —Sa coincide only in this 
limit. The evolution equation of Qsm reads [46, 52] 



H \ a H a' I a z 



Qsm + 3-ff Qsm + 

On large scales the adiabatic evolution is simpler and governed by the first integral 



sm =-2(f)',)s)' -■2f)'[^- + ^)6,. (75, 



Q'sM+[jj-^r)QsM-2e'8s^Q, (76) 

which simply expresses that Se, written in terms of Qsm and Ss, is negligible on large scales. One can easily check 
that the large-scale limit of (75) follows from the first integral (76). 

Let us now consider the evolution equation for £ a . The spatial components of £ a , at linear order, are given by 
cq. (18) as for the fluid case: Q = di(. Similarly, one can define 1Z by IZi = dilZ. Note that 1Z coincides with the 
familiar comoving curvature perturbation only in the large scale limit, because of the difference between a and ip, as 
already discussed for (. Eq. (66) implies a simple relation between ( and 1Z, 

C + K=~6e, (77) 
P 

which shows that £ and —1Z coincide on large scales. 

The evolution equation for ( on large scales follows from eq. (67) which, after linearization, yields 

C « -™e'5s. (78) 
a 

One thus recovers the familiar linear result [46] that the entropy perturbation is sourcing the evolution of the uniform 
density curvature perturbation on large scales. 



2. Second- order perturbations 

As for the perfect fluid case in section II B 2, we now expand the equations governing a a , s a and ( a at second order 
in the perturbations. For simplicity, we will restrict ourselves to large scales and we will thus start from the equations 
expanded in spatial gradients discussed in Sec. IIIB. 

The second-order evolution equations for the perturbations of a single scalar field, in a coordinate based approach, 
has been considered in several references (see for instance [11, 12, 53, 54]). The multi-field case has been first studied 
by Malik in [55, 56] and, using the separate universe approach, in [57, 58]. Here, we obtain directly the second-order 
evolution equations in terms of the adiabatic and entropic perturbations. 
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By expanding Eq. (47) and Eq. (48) up to second order, one can write [25] 

5a (2) =d t 5a (2) + ^Sad t Ss - 



5a, 



with 



§s W =diSs W + »JL di§s ', 



La 



3 (2)=_^^(2) 

and where we have defined the spatial vector 



5a 



a' a' a' 1 



+ 77^ > 



Vi = -(5sdi5s' -Ss'diSs). 



(79) 
(80) 

(81) 
(82) 

(83) 



The definition of Ss^ is chosen such that it is gauge invariant on large scales. Since the adiabatic component a a does 
not vanish at zeroth order, 5a ' is not a gauge invariant variable and can be chosen for convenience. Its definition 
here is such that the momentum perturbation almost vanishes (up to Vi) when 5a^ = 0. 

The presence of the spatial vector Vi on the right-hand side of Eq. (79) is due to the fact that a a defined in 
Eq. (47) is not hypersurface orthogonal. Indeed, one can check using its definition that cr[ a Vfc(T c ] = a a^ a SbS c j . As 
the momentum 5qi is proportional to 5a^~, Eq. (58), this implies that if Vi does not vanish one cannot define at 
second-order a comoving gauge, i.e., such that Sq^ = and 5q^ = 0, in contrast with the linear theory or the 
single- field case [48]. However, the evolution equation for Vi on large scales reads 



V' + 3HVi « 0, 



(84) 



which implies that in an expanding universe Vi decays as a 3 and rapidly become negligible even if it is nonzero 
initially. Consequently, in the following we will ignore Vi on large scales. 

Since 5a {2) is not gauge- invariant, it is useful to consider the Sasaki-Mukhanov variable at second order Qg^, which 
on large scales coincides with the scalar quantity which can be extracted by expanding at second order Q a defined 
in Eq. (62), 



Q {2) = d t Q^ + ^Q' + 6 -Qd t 5s - its. 



The Sasaki-Mukhanov can be found by replacing Eq. (26) in the definition of Q^ 2 \ which yields 



lfa'Y 



2\H 



V> - 9' 5s 



(85) 



(86) 



Restricted to a single scalar field, this definition coincides with the one given in [55]. By using the energy and 

momentum constraint equations at second order it is possible to derive a first integral for [25], extending at 

second order Eq. (76). 

Here we will simply derive the large-scale evolution equation for £ at second order, which is enough to describe the 
evolution of the adiabatic mode. Expanding up to second order Eq. (67) yields 



c 



(2)/ 



H 

^2 



29'a'5s {2) - (V ss + AO 12 ) 5s 2 + ?f 5s6s' 



It is also useful to express our results in terms of lZ a . The spatial components of lZ a can be decomposed as 

a' a' z 



with 



= - Sa W + -5a^ + - 
a' a' 



1/iA' 



2 \a' 



Sa + O'—Ss 
a 1 



(87) 



(88) 



(89) 
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The last term in Eq. (88) comes from the fact that, in contrast to ( a , lZ a is defined in terms of the spatial momentum 
which cannot be expressed in general as a pure gradient. However, in an expanding universe this term can be neglected 
and IZS 2 ^ coincides with the second-order comoving curvature perturbation defined in [54, 59]. 

It is easy to derive a first-order (in time) evolution equation for VS 2 ^ by noting that and TZ^ are related on 
large scales. Indeed, by using Eqs. (24) and (88), neglecting Se^ and Vj, and making use of (p) + IZ^ « 0, one 
can show that the spatial component of Eq. (66) yields £( 2 ) + VS 2 ^ « [54]. From this relation and the second-order 
evolution equation of £, Eq. (87), one can find a large-scale evolution equation for 1Z at second order, 



29'a'8s (2) - (V ss + AO' 2 ) Ss 2 + ^-SsSs' 



(90) 



The second-order perturbation can be related on large scales to by combining Eqs. (86) and (89) and using 
Eq. (26). One obtains 



Q 



(2) 



1 fH 



> m ~ ^7 (Q'sm - & 6s) Q SM - — f — ) Qsm 



(91) 



(2) 

which can be used to show that Eq. (90) can be rewritten as a first integral for Q^m- 

Let us now discuss the second order evolution of the entropy perturbation Ss. On large scales this is obtained by 
simply expanding the spatial components of Eq. (65) up to second order, which gives 

Ss^" + 3H5s™' + (V ss + 39' 2 ) 5s™ « -^5s' 2 1 (V + 9'^ - ^hA SsSs' (W sss - 5^V SS 9^) Ss 2 . 

v ' a' a' \ a' 2 / \2 a' & ) 

(92) 



As in the linear theory, the entropy perturbation evolves independently of the adiabatic component on large scales. 



3. Third-order perturbations 

It is straightforward to extend the same procedure at third order, as discussed in detail in [43]. The evolution 
equation for £( 3 ) on large scales can be obtained by expanding up to third order Eq. (67), which yields an equation 
of the form 

^Sf[Ss,Ss',Ss^,Ss^'], (93) 

where the source term on the right-hand side contains terms cubic in first order perturbations (Ss or Ss') and terms 
that are products of second order perturbations (Ss^ or Ss^') with first order perturbations (Ss or Ss'). The explicit 
definition of 8s^ used in this expansion is 

Ss^ ee e sI S<f>^ - ^ (Ss' + 9'Sa) %Ss^' - ^(Ss" - ^Ss' + 9 2 Ss) ^SsSs'Sa, (94) 

u' 'a' 2c z a ocr \a J 2a 1 

which can be determined by an appropriate expansion of s a at third order, similar to the second order expansion (80). 

To complete this equation, one needs to derive the second order (in time) equation of motion at third order on large 
scales. This is done by expanding Eq. (65) up to third order, which yields an equation of the form 

«Js< 3 >" + 3HSs^' + (% s + 39 2 )Ss^ w [6s, 5s' ,6s™ ,Ss^'}. (95) 

The explicit expressions of the source terms and S s , which have been derived explicitly in [43], are given in the 
Appendix. 



D. An application: super-Hubble generation of non-Gaussian perturbations 



An important application of the above formalism is the computation of the non-Gaussianities of the primordial 
curvature perturbation. The computation of £ and Ss, by integration of their respective evolution equations up to 
the required order, can be seen as an alternative to the SN formalism, described in the article by Wands in this 



13 



issue [60] . It is especially useful when the background evolution is complicated so that one cannot express analytically 
the number of e-folds as a function of the initial scalar field amplitudes. 
Formally, one can write the solutions of the evolution equations as 

C = C* + T^Ss* + if <5s 2 + if 6s 3 + ... , (96) 
5s = T^Sa. + T^Ssl + T^Ssl + . . . , (97) 

where the if and 1'j™' 1 correspond to transfer functions, for £ and 5s, respectively, at the n-th order and £*, 5s* are 
their initial conditions at Hubble exit. In particular, we are interested in the values of £ and 5s after inflation, at the 
onset of the radiation dominated era. 

Neglecting slow-roll corrections, the field perturbations Qsm* and 5s* at Hubble crossing can be treated as Gaussian, 
with respective power spectra PQ,(k) and P$ Sit (k) (see Appendix of [61] for an explicit calculation of their 3-point 
functions). Note that there exist models with non-standard kinetic terms where the perturbations just after horizon 
crossing are non-Gaussian (see [7] and the article of K. Koyama in this issue [62] for a general discussion on multi-field 
models of this type) and where the entropy perturbations can play a crucial role [63]. However, we will not consider 
these models here. 

The perturbations Qsm* and 6s* are also independent random fields, at least at leading order in slow-roll [64, 65]. 
Moreover, although the relation between £ and Qsm is in principle nonlinear (as illustrated for example by Eq. (91) 
for 1Z which coincides with £ on large scales) , at horizon crossing the linear relation dominates in slow-roll models so 
that C* ^ -(H/a')Q SM * [61]. 



1. Curvature perturbations 

Using Eqs. (96) and (97) one finds that the final 2-point, 3-point and 4-point functions of £, which are in principle 
observable, are given respectively by 

(CfeCfe,) = (2nf6(k + k')P c (k), P c (k) = (2e*M 2 )- 1 P Q ,(k) + (T^) 2 P Ss ,(k), (98) 

(Cfe, k 2 Cfe 3 ) - 2(T C (1) ) 2 T C (2) [P Ss , (k^Pss, (k 2 )+2 cyclic] , (99) 
(C^CkAs^J - 4(T C (1) ) 2 (T C (2) ) 2 [P*,.(fci)P 4a .(fc 2 )P* a .(fci3) + 11 Perms] 

+6(T C (1) ) 3 T C (3) [PssAki)P5sAk2)PssAk 3 ) + 3 cyclic] , (100) 

where we have introduced the slow-roll parameter e = inGa' 2 / H 2 and Mp = (87rG)~ 1//2 . 

If the power spectrum is dominated by its entropy contribution, i.e. (T^) 2 Pg St ^> (2e*Mp)^ 1 Pq^ , then one can 
easily relate the non-Gaussianities generated by the isocurvature field to the local nonlinear parameters /nl, tnl and 
5nl, defined by [57, 59, 66, 67] 

(CkAAJ = !|/nl [P C (fci)P C (fc 2 ) + 2 cyclic] , (101) 

54 

(CkAAAJ = t nl + 11 Perms] + - 5NL [P c (/c 1 )P c (fc 2 )P c (fc 3 ) + 3 cyclic] . (102) 

From Eqs. (99) and (100) one finds 



Tf (if ) 2 25^ 3) 

3 (T (D )2 ' TNL - 4 (T (D )4 ' ffNL " 9 (T (i) )3 



/NL — 7T 



t nl = 4^^, ffNL = -T^iw (103) 

c 



These expressions can be applied to compute the non-Gaussianities generated in multifield models. For instance, 
in the case of two fields, for separable potentials - potentials that can be written as the sum or product of two 
functions dependent only on one of the fields, the transfer functions can be computed analytically either from the 
evolution equations for adiabatic and isocurvature perturbations derived in this section or by using the 5 TV-formalism 
[58, 68, 69]. It has been shown that the non-Gaussianity generated after inflation are slow-roll suppressed in these 
cases (see also the review by Byrnes and Choi [70] ) . The non- linear evolution of the isocurvature field during inflation 
can generate a large non- vanishing 3 and 4-point functions [71, 72], as reviewed by Bernardeau in this issue [73]. From 
Eq. (103), the conditions to obtain large nonlinear parameters are if > (Tf) 2 and if > (if) 3 . 
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Eqs. (103), together with the evolution equations for adiabatic and isocurvature perturbations derived in this section, 
have also been applied to compute the non-Gaussianities generated during the new Ekpyrotic scenario [43, 74-76] 
with results that agree with those obtained using the SN- formalism [43, 77]. 

A similar approach (although keeping the original scalar fields instead of a decomposition into adiabatic and entropic 
modes) has been used to compute the evolution of £ up to second order, and thus obtain the non-Gaussianities, in 
the context of hybrid inflation [78-81]. In these works the curvature perturbation at second order was expressed as a 
time integral of an expression quadractic in the linear perturbations of the tachyonic field. The same procedure can 
also been used to study the preheating phase at the end of hybrid inflation [80, 81]. 

2. Isocurvature perturbations 

It is possible, although not necessary, that the entropy perturbations generated during inflation survive after 
inflation, in the form of traditional entropy or isocurvature perturbations, i.e. fluctuations of the relative number 
of particles between two species (for instance, the number of photons per cold dark matter particle). In this case, 
the accessible cosmological data would result from a combination of curvature and isocurvature perturbations, which 
could be correlated [82]. So far, there has been no detection of an isocurvature perturbation in the CMB data and the 
constraints on the power spectra are already tight [83] . One could also envisage that an isocurvature component might 
be detectable through its own higher-order correlations or its higher-order correlations with the adiabatic component, 
as discussed in [61, 84-87]. Such a possibility can arise for instance in the mixed curvaton and inflaton scenario 
introduced in [88]. 

In the radiation dominated era the adiabatic perturbation coincides with the perturbation in the radiation fluid 
which, using Eq. (12) with w = 1/3, can be written as 

C a = d a (, C^Sa + ]ln(^). (104) 



4 \P, 

The isocurvature perturbation between the cold dark matter and radiation, say, can be defined as 

Sa = 3(C - C) = d a S, S = In (|) - \ In (|) . (105) 

Formally, one can write the relation between the entropy perturbation during inflation and the final isocurvature 
perturbation in the form 

S = Tf ) 5.^+Tf ) 8sl+Tf ) 8sl + ... , (106) 

where the transfer functions Tg 1 ^ depend on the details of the model. Analytically solvable examples of Tg 2 ^ are given 
for instance in [61]. 

The observational upper limit on the fraction of isocurvature fluctuations allowed by data (see [83]) implies that 
the linear transfer coefficient Tg 1 ^ must be very small. Interestingly, however, t| 2 ^ and Tg 3 ^ could still be large and 
the non-adiabatic perturbation detectable through its effect on higher-order correlation functions. Indeed, neglecting 
Tg~\ one finds that the non- vanishing 3 and 4-point correlation functions are given by 

(Ck^A) * 2(T^) 2 Tg 2) P Ss , (k^Pss, (fe), (107) 
(CkSA^kJ - 4(T c (1) ) 2 (Tf ) 2 P 5s , (h)P Ss , (k 2 ) [P S „ (fc 13 ) + P Ss , (fc 23 )] , (108) 

(CkAAs S kJ * 4(r c (1) ) 2 r c (2) r| 2) {p, St (h)p Ss , (fc 2 ) [p s „ (k 13 ) + p Ss , (fe 3 )] + 2 c yc iic(i, 2, 3)} 

+ 6(T^ } fTg 3} P Ss , (h)P Ss , (k 2 )P Sa . (A3). 

It would thus be interesting to obtain some observational constraints on the isocurvature non-Gaussianities from the 
future data, such as those collected by the Planck satellite. 
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Appendix: source terms at third order 

The source terms that appear on the right hand side of (93) and (95) are respectively given by [43] 

29'a'Ss^ - 2{V SS + 46' 2 )6s5 s W + ^-(Ss^SsJ 

1 



and 



S^ 3) [Ss,8s',Ss^,8s^'} = -^ 



+ 



3 & 



2 (-„„ ^,V„ 3 



12 



V 
a 

8' 3 9' 
8— +4— V 



,ss 

a' a' 3 



a' 
V. 



H& {SsSsWy - v, 



s.-S — 100' 



1 8 1_ SsSs^ 



vi 



3^ + 3H^r - - 6^ Ss' 2 6s 
a 6 a 1 a 1 ] 



a + 7—7^ + 3i7^f - 14i?^ 5s' 5s 2 



a ' 

-/. 
(7 ' 



(T ' 



(7 

,"'4 ' 



5' 2 



~V;„ a + 24f + 21^j^ + 21 u — ) fc 3 



(110) 



(111) 
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